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ABSTRACT 


A study of the vapor-liquid phase equilibrium 
for the system Nitrogen-Methane-n-Hexane was made at 
temperatures of 40° and 100°F and pressures of 500, 

1,000 and 2,000 psia. 

A variable volume, windowed equilibrium cell, 
with mercury as the confining liquid was used in this 
work. The experimental and analytical equipment and 
techniques used are described. 

The analyses of the equilibrium phases at the 
different conditions of temperature and pressure are 
presented. The equilibrium ratios were calculated using 
the experimental data and were compared with those cal- 
culated using the Winn Nomograph and the NGPSA Engineering 
Data Book. It was found that these methods do not give a 
good prediction of the K-factors, probably because of the 
azfficulty in making a reliable calculation Of the con- 
vergence peesetre at the temperatures investigated for 
this system. 

The Chueh-Prausnitz method was used to predict the 
critical loci of the Nitrogen-n-Hexane and Methane-n- 
Hexane systems. For this last system, no critical experi- 
mental data were used in the prediction, but the calculated 
locus compared very closely with that experimentally ob- 
tained. 

The Chueh-Prausnitz method was used to calculate the 


bubble point temperature and vapor composition from a 
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knowledge of the pressure and liquid composition for 
the Nitrogen-Methane-n-Hexane, Nitrogen-n-Hexane and 
Methane-n-Hexane systems. When compared with the 
experimental results, it was found that the prediction 
of the bubble point temperature and vapor composition 
were not accurate. The binary data needed for the 
calculation had to be obtained by extrapolation, and 
it is believed that the available data were insufficient. 
It was concluded that before reliable calculations 
could be made using this method, experimental informa- 
tion on the binary pairs must be available over the 


entire temperature range. 
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fap DL RODUCT ION 


The vapor-liquid equilibrium behavior of hydro- 
carbon and hydrocarbon-non-hydrocarbon systems is of 
great practical and theoretical importance. 

The properties of coexisting phases are needed 
for the proper design and operation of industrial equip- 
ment encountered in the chemical-petroleum processing 
field such as distillation columns, absorbers and evapo- 
FakroOrs. 

Experimental and theoretical methods are now used 
to describe the properties of the equilibrium phases in 
hydrocarbons and related systems by specifying the equi- 
librium ratios of the components. The theoretical methods 
used are either thermodynamically or empirically developed. 
Assumptions made during the formulation of those methods 
Limite thesrange.of stheur application. 

The vapor-liquid equilibrium of many hydrocarbon- 
non-hydrocarbon systems have been investigated, but there 
are still many systems of practical and theoretical an— 
terests that have not yet been studied. The prediction of 
Phesequidibrium,ratios, ter these systems ghas, vecentily been 
facibitatedsand amproved byemeans Ofenewly developed me- 
thods such. as.the, Chueh-Prausnitz method, This method 
has the advantage that the effect of the composition on 
the activity coefficient of the liquid is calculated 


separately from the effect of the pressure and that only 
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properties of pure components and binary pairs are 
needed for the prediction of ternary and multicomponent 
phase behavior. 

The fact that nitrogen is present in appreciable 
amounts in petroleum reservoirs, made the study of the 
phase behavior of nitrogen with hydrocarbon systems very 
important for the production and processing of the petro- 
leum fluids. As long as thirty years ago, Boomer and 
Johnson, at the University of Alberta, realized the 
importance of the non hydrocarbon-hydrocarbon systems. 
They Originated the study of the effect of nitrogen on 
hydrocarbon systems by investigating the vapor-liquid 
equilibrium of the following systems: Nitrogen-Methane- 
n-Pentane, Nitrogen-Methane-(ntiso) -Pentane, Nitrogen- 
Methane-n-Hexane, Nitrogen-Methane-(n+iso)Hexane, 
Nitrogen-Methane-n-Heptane. They used a gas composed 
of 5.6% nitrogen and 94.4% methane, and studied the 
Nitrogen-Methane-n-Hexane system in the range of tempera- 
ture from 25°C to G5 Cc, and the pressure from 35 to 230 
atm. 

It was only recently, that Poston and McKetta, 
studied the same system at temperatures from 100 to 
340 F and pressures from 500 to 5,000 psia. 

For this work, it was decided to continue the study 
of the Nitrogen-Methane-n-Hexane system at a temperature 


of 40°F and at pressures up to 2,000 psia. 
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The objectives of this study were 


1. To extend the range of temperature over which 
K-factors were available for the Nitrogen-Methane-n- 
Hexane system. 

2. TO compare the K-factors obtained by other methods 
like the NGPSA and the Winn Nomograph, at the high con- 
vergence pressure caused by the presence of nitrogen. 
3. To try to reach a conclusion about the best method 
available for the prediction of the K-factors and im- 


prove the prediction if possible. 
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Ti LAH HORY 
The Phase Rule 


The phase rule given in 1876 by J. Willard Gibbs 
is still frequently used in phase behavior studies. It 
makes it possible to calculate the number of variables 
that have to be specified in order to completely define 
a system at equilibrium. It can be mathematically ex- 


pressed as follows: 
Hos Cant Pite2 Gl) 


where F = the number of degrees of freedom of the system 
C = the number of components of the system 
P = the number of phases in equilibrium in the 
system. 
Therefore in a three component system, when a state 
of equilibrium exists between two phases, three independent 
variables have to be specified in order to completely 


describe the system. 
Criterion of Equilibrium 


The necessary conditions for equilibrium to exist 
between any mixture of components in the vapor and liquid 


phases are 


pv = oh (2) 
by eae (3) 
ta feert (4) 
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where T', T are the temperatures of the vapor and 
liquid phases 
iy G2 are the total pressures of the vapor and 
liquid phases 
: ie are the fugacities of component i in the 


vapor and liquid phases. 


Equilibrium Ratio 


The equilibrium phase distribution ratio or K-factor 
is the most common equilibrium quantity used in process 
calculations for coexisting vapor-liquid mixtures. It is 


arbitrarily defined as 


e al 
ie We (5) 
al. 


where Y; and X. are the mole fractions of component i in 
the vapor and liquid phases. 

Several texts (45, 46, 22, 16) have dealt with the 
thermodynamics of phase equilibria at high pressures. 
These books have included detailed discussions on the 
different methods available for predicting the equilibrium 
ratio for hydrocarbon and other components of interest to 
the petroleum industry. No method is completely general; 
every method involves certain assumptions which make it 


applicable over a limited range. 
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Ideal Equilibrium Ratios 


At low pressure, the fugacity can be replaced by 
the pressure and the criteria of equilibrium could be 


taken as 


= )p. (6) 


where Das and Bae are the partial pressures exerted by 
component i in the vapor and liquid phases. 

Assuming that the coexisting vapor and liquid 
phases are both ideal solutions, that the vapor phase 
behaves as an ideal gas, and that the vapor pressure of 
the liquid is a function of temperature and is not affect- 
ed by the pressure; then a combination of Raoult's and 


Dalton's law could be used to calculate: the equilibrium 


ratio 
Se Pp. Raoult's law (7) 
i ia 41 
Sey : 
P; ag ee. Dalton's law (8) 
Ys Po? 
At equilibrium K, = xa Ferre (9°) 


where PS is the vapor pressure of component i at the 
temperature of the system and P is the total pressure of 
the system. 

Generally speaking, K-factors predicted by this 
method are restricted to low pressures, far below the 
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Ideal Solution Equilibrium Ratios 


To eliminate the assumption made in Dalton's law 
of ideality in the vapor phase, the fugacity could then 
be used instead of pressures. With the assumptions 
that the liquid and vapor phases are both ideal solutions, 


the Lewis-Randall fugacity rule could be applied 


Le = xh £0 (10) 


ey ee Cyl es (11) 


where fe is the fugacity of component i as a liquid at 
the pressure and temperature of the equilibrium 


fe 
di 


is the fugacity of component i as a vapor at 
the pressure and temperature of the equilibrium. 

Many methods have been reported to calculate the 
ideal solution equilibrium ratio. The ideal solution 
equilibrium ratio is only a function of the temperature 
and pressure of the system and does not depend on the 
composition. 

The fugacity of either the liquids or vapors are 
sometimes obtained by extrapolation of the fugacity 
pressure relation to a hypothetical state. The ideal 
equilibrium ratio could also be calculated using an equa- 
tion of state that represents the PVT behavior of the 


components in both the liquid and the vapor phases of the 


system. It has been found that the equilibrium ratios 
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obtained from the ideal solution theory do not hold as the 


Critical conditions are approached. 


Actual Equilibrium Ratios 


The substitution of the activity for mole fraction 
in the equations derived for an ideal system, made those 


equations valid for real systems 
(12) 


si yee yore (13) 


where i and vi" are the activity coefficient of compo- 
nent i in the liquid and vapor phases respectively at the 
system temperature and pressure. 

Therefore, all the deviations from ideal behavior 
are corrected by the introduction of the activity coeffi- 
cient which changes with temperature, pressure and com- 


positron. At-equrlibrium 


y L L 3 L 
a i Z 2 (14) 


There are no assumptions involved in formulating this 
equation but the difficulty is in finding a general 


method for calculating the activity coefficients. 


Convergence Pressure 
From the experimental results obtained for differ- 


ent systems, it was observed that the K-factors converge 
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to unity at some pressure for each temperature chosen. 
This pressure is called the convergence pressure. 

The convergence pressure of a system is thus de- 
Ginedtas that pressurerterpuwhich ~EhesK-factors of all 
components in the system converge to unity at the system 
temperature, 

Hanson and Brown (23) showed that for paraffin 
hydrocarbon systems that for a given component, the 
equilibrium ratios are the same for the same convergence 
pressure, even when the composition of the system is 
different. 

The convergence pressure has been widely used as 
a composition parameter, to relate K-factors to the com- 
position of the hydrocarbon systems. 

When applying the phase rule to a two phase system, 
the number of degrees of freedom is equal to the number 
of components. If the convergence pressure is expressed 
as a function of temperature and N-2 composition parameters, 
the system wills be completely described by fixing the pres- 


sure and the convergence pressure, then 


1 SDB, Ee (15) 


For a binary system, the convergence pressure 
depends on the temperature alone, and is equal to the 
critical pressure of the system that has the equilibrium 


temperature for its critical temperature. 
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Several methods have been suggested for the pre- 
diction of the convergence pressure of multicomponent 
systems, and for the calculation of the equilibrium ratio 
using the convergence pressure concept. 

Hadden (20) showed that the convergence pressure 
is a function of liquid phase composition excluding the 
concentration of the lightest component. He suggested 
the use Of a composition parameter based on mass 
fraction to obtain the convergence pressure. For a tern- 


ary system, the composition parameter would be 
m 
M, = ————-. (16) 


where the components are numbered in the order of de- 


Creasing, volatility. 


m, = mass fraction of component 2 in the liquid 
phase 

m, = mass fraction of component 3 in the liquid 
phase 


Using the composition parameter suggested by Hadden, 
Winn (78) developed an expression for the convergence 


Pressure of a ternary mixture as follows 
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P. and P, are the convergence pressures for the binaries 
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composed of components (1, 2) and (1,3) respectively at 
the system temperature. 

If the system temperature is greater than the 
critical temperature of component 2, then P, is taken 
equal to.the critical pressure for ;icomponent 2; 

PL is the convergence pressure for the ternary 


mixture at the specified values of M, and at the tempera- 


2 
ture of the system. 

He presented a mnomograph, ito, obtain |the kK-factors 
when the convergence pressure was 5,000 psia. K-factors 
could be estimated when the convergence pressure was other 
than 5,000 psia, using a grid pressure evaluated from 
another chart. 

Rzasa et al (56) devised a correlation to estimate the 
convergence pressure of complex hydrocarbon systems. They 
based their correlation on the relationship of the critical 
pressure of the system and the product of the molecular 
weight and specific gravity of the heptanes plus fractions. 
The obtained convergence pressures are then used to get 
the equilibriumsrataos. 

Hadden and Grayson (21) extended the Winn nomograph 
to obtain the equilibrium ratios at temperatures from -260°F 
to 800°F and at pressures up to 10,000 psia. These nomo- 
graphs apply to light hydrocarbons and petroleum fractions 
and to certain non-hydrocarbon systems. Nitrogen is inclu- 


ded in thosenomographs. The K-factors are related to 


i'8 Aire - ye 4 = SUE hee 
- et eee 
a Tes seme Lee 


Se ; a 
ne vi) bas (Sf) etnenognos. 20.58 OK 


: ee 
: tata; | | .oiussi9qmes moses 


ed neit s9t6970 et sivtsxegtet meteye sds tr 
moras et of aeit ,S tnamogmoo to sirvssisqmes Iso£ 


«§ snpmogmos x03 Ssivesoid {s0isixo sid oF te 


yxsazet oft 10lt siveaotq sonaptevron eft at at 


~si9sqmes edt 35 Bris ol Ro zeulsv helifosge sft +5 
sneteye oft Xo exud 

erotost-aA edt otstdo ot dqsipomon s betasesitd 5H 
atotost-N .sia@q 000,c esw e1wessig soneprsva09 sat as 


‘tedto e6w emeesiq soneprovaoco et nedw betsmisas sd bivoo 
- , a 
mox? bestsulsvs siueesig bixp s pnieu ,Siaq 000, @ asda 


_txsdo xeds 5 
~t 7 r 7 
' - 


eft stsmices o+ noitsisizr0o> 5s bsatveb (ec) Is js sessaA a 
yedT .emetaye nmodxssoubyd xelqmoo to sruaesiq Ssonept 
Isvttinzo edt to dindenoitsiex els no noitsisitoo tient be 


sxeluoolom sds to touborq efit bas meteye sdt to otees? ce 
R i 


-anoisosx? aulq esnstged sist io ytivsip oLiiLosge bas sapien 
jsp ot been aesdt 916 asuveestg SonsprevnoD bentsado eft 

peal ,2oftsx ing ach edad 
riqsxpomon aniwW sft bebaetxs (iS) noeysx9 bas sopbes 


a°0as- moxt SeSOSETaES 35 eckss “initia te eds mis wall ) O3 
peta 
~omon eeenT .ateq 000,0L ot qu eerveesig te brs 3008 « 


“eaoksoss? malozseg bas enodissorbhyd tiytt o3 vias edger 
oer wl aduduied : cabal dodtsvoxbyd-non nistis of 
Se seth 
vF i od betsl nt os fenton 2 Sa ho BX ont A 
yee xh “a c Ti a eon ons 7 ad 


ae 


- 12 - 


component temperature, pressure and convergence pressure. 

Myers and Lenoir (35) presented a nomograph to 
evaluate the(|K-factorsausing -thetvalve+o& Keatel0'psia 
and a convergence pressure of 5,000 psia. 

The values of the equilibrium ratios are increas- 
ingly affected by the convergence pressure when the 
operating pressure approaches the convergence pressure. 
Therefore, the convergence pressure should be accurately 
calculated to get a proper estimate of the equilibrium 
ratio. 

The convergence pressure methods are the most 


widely used empirical methods to obtain the K-factor. 


Chao-Seader 


Chao and Seader (11) developed a correlation of vapor 
liquid equilibria in hydrocarbon mixtures.’ They calculated 
the equilibrium ratio for the components in the equilibrium 


mixture using a three parameter equation 
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where se is the activity coefficient of component i 
in the liquid solution 
ee is the fugacity coefficient of pure liquid 
component i at the conditions of the system 
o.° is the fugacity coefficient of the component i 


in the vapor mixture. 


The Redlich-Kwong equation of state was used for 
the calculation of the fugacity coefficient of component 
iin the vapor mixture. To apply the equation, two con- 
stants are required for every component, these are directly 
related to the critical temperature and pressure. 

Chao and Seader calculated the activity coefficient 
in the liquid with the assumption that the liquid hydrocar- 
bon solutions are regular solutions. Regular solutions have 
their excess entropy equal to zero, any non-ideal behavior 
is due to the heat of solution. They used the equation 


given by Hildebrand 


V. (6,5 - ype 
Ral 


where v. is the liquid molar volume of component i 
i 


6 is the solubility parameter of component i 


The solubility parameter is equal to the square 


root of an energy density i.e. 
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where AE. is the energy of vaporization 


(21) 


The solubility parameter and molal volume could be 
taken as functions of temperature. However, it was 
found that they could be considered as characteristic 
constants of a component without appreciable loss of 
accuracy. 

The liquid fugacity coefficient of pure compo- 
nents have been correlated with a modified form of the 


principle of corresponding states given by Pitzer. 


Accordingly 
Hognesh= (toy joe oneToeere (22) 
ee) and have been fitted with approximating func- 


tions and depend only on reduced temperature and reduced 
pressure. 
w is the acentric factor, it gives. a measure of 

the deviation of the behavior of usual substances from 
Phatweo beanvuidealizedasimplé,fluidveamlters sancharacteristic 
constant for each component. 

The Chao-Seader correlation can give high accuracy 
in the range specified below: 
1. For hydrocarbon (except methane), a reduced temperature 


range from 0.5 to 1.3 based on the pure component Cplcican 
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SUS ae 
temperature and pressures up to 2,000 psia but not to 
exceed about 0.8 of the critical pressure of the system. 
2. For systems containing light gases like hydrogen and 
methane 
a) The reduced temperature should not exceed 0.93 
b) The temperature lies between -100 F and 500;-F 
c) Pressures up to 8,000 psia 
d) Other dissolved gases should be below 20 mole 


percent in the liquid. 


Lenoir and Koppany (33) showed that for binary mixtures 
of methane in hexane and paraffin of higher molecular 
weight, the Chao-Seader method predicts the equilibrium 
ratio of the heavier paraffin very poorly. 

Prausnitz (44) mentioned some limitations of the 
Chao-Seader correlation and stated that this correlation 


is unreliable at conditions approaching the critical. 


Chueh-Prausnitz 


Prausnitz and Chueh (45) presented equations for 
reducing binary vapor-liquid equilibrium data to thermo-— 
dynamically significant binary parameters; upon generaliz- 
ing these equations to systems containing any number of 
components, they were then able to predict the phase 
behavior of multicomponent mixtures. Their method is 
restricted to non-polar or slightly polar fluids as those 


encountered in the petroleum industries. 
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For the vapor-liquid equilibrium of an N-component 
system, there are N degrees of freedom. 


For a system in equilibrium 


Wiese Sisus ple 2, 3 ek EOUN 


Prausnitz and Chueh made use of the following equations, 


to be able to solve the N simultaneous equations given 


by (4) 
ee =e Ye (23) 
ce ee Te aa 
ee - 1, exp i. a ap (25) 
P 


where ds is the fugacity coefficient of component 12 in 

the gas mixture 

P is the total pressure of the system 

£2) por the solvents, it is the fugacity of pure 
liquid i at the temperature of the solution and 
at the constant reference pressure ae |, 
For the solutes, it is Henry's constant for the 
solute i in the solute-free solvent mixture at 


system temperature and at the constant reference 


c 
pressure P-. 
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(ele a oe 
i is the activity coefficient for component i 
a Gee 
(Bied. ie vs | 
ia 1s the activity coefficient for component i 
at the reference pressure oo 
ba 
1B Standard State Pressure 
Vv Partial molal volume of component i in the 


Piquid mixture. 


The simultaneous solution of the equations given by (4) 
is most conveniently performed by an iterative procedure. 


The solution must satisfy the two stoichiometric relations 


and 


A thermodynamic description of high pressure vapor- 
liquid equilibrium must describe how the fugacity of each 
component, in each phase, depends on the temperature, 
pressure and composition. 

In the vapor phase, this dependence is given by the 
fugacity scoefficient,, and .in the liquid phase; this iWde= 
pendence is determined essentially by the effect of tem- 
perature, pressure and composition on the activity coeffi- 
Certs. 

It is necessary then to calculate the fugacity co- 


efficients, the activity coefficients, the reference 
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fugacities and the partial molar volumes. 

ds was calculated as a function of (P, T, Y;) 
using a revised Redlich-Kwong equation of state (to 
take into account the deviations from the geometric 
mean approximation) with new mixing rules. 


E 


Pp” the standard State pressure was taken for 


convenience equal to 0 


It is an experimentally accessible quan- 
ELty 

Ve was calculated in terms of T, X, and the satu- 
rated liquid molar volume of the mixture. The Redlich- 
Kwong equations of state was applied to fit the PVT data 
of saturated liquid for each pure component. Modified 
mixing rules were used to apply the equation to the liquid 
mixture using corresponding states correlation. 

The symmetric convention of normalization was used 
for the activity coefficients of solvents, and the un-— 
symmetric normalization was used for solutes. 


ie 
oo ) was calculated as a function of T and xX, using 


a two parameter equation obtained from a dilated van Laar 
model. To obtain the constants of the equation, self- 
interaction and dilation constants are needed, those are 
obtained from binary data. The van Laar's equation was 
modified to allow for the rapid change of the liquid molar 
volume which occurs in the critical region. 


The method can be used to calculate multicomponent 
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vapor-liquid equilibria for non-polar or slightly polar 
systems, but in addition to the pure component properties 
some binary data for all possible binary ypairs_in,.the 
mixture are needed. 

The method has not yet been extensively tested in 
practice but it is considered to be one of the best pro- 
posals to date for predicting vapor-liquid equilibrium 


at high pressure. 
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III RELATED SYSTEMS INVESTIGATED 


The properties of pure nitrogen and pure methane 
have been extensively studied and are well established. 
Less work has been done with n-hexane. The properties 
of n-hexane at atmospheric pressure have been studied in 
some detail and were summarized by Rossini (55). The 
earlier work of Shepard (67) served to establish the 
index of refraction and the specific weight at atmos- 
pheric pressure which was also evaluated by Wibaut (76, 
77). Kelso and Felsing (27, 28) studied the volume- 
temperature relations of n-hexane at pressures up to 
4,500 psia, and Eduljee and co-workers (17) determined 
the volumetric behavior of this compound at pressures up 
to 75,000 psia. Stewart, Sage and Lacey (71) investi- 
gated the volumetric behavior of n-hexane at pressures 
up to 10,000 psia and temperatures between 100° and 460°F. 
Leverett (34) determined the PVT properties of n-hexane 


at 200-800 mm Hg and 40°-200°C. 


Related Nitrogen Binary Systems 


The vapor-liquid equilibrium of the nitrogen-methane 
system has been studied by several investigators (9, 13, 74, 
36, 18). McTaggart and Edwards (36) studied this system 
at atmospheric pressure and their results show large devia- 
tions from ideal solution behavior. Torochnesnikov and 


Levins (77) determined six isotherms for the nitrogen- 
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methane system between -298° and -200°F, while Cines et al 
(13) made their investigation in the temperature range of 
“200 F and -150°r. "Bloomer ana Parent (7) studying the 
same system, covered the range of pressure from atmospher- 
ic through the critical region. The lowest temperature 
they studied was -295°F, and the highest was the critical 
temperature of methane -116.7°F. They reported some ex- 
perimental values for the densities of the saturated 
liquid and saturated vapor phases of pure methane, and 
were able to draw the critical locus curve. The agree- 
ment was in general excellent between the data reported 

by Bloomer and Parent, and Cines et al. The greatest 
difference was in the critical region. The data of Bloomer 
and Parent, and Cines et al, do not agree with the earlier 
data reported by Torocheshinkov and Levins. 

Akers, Attwell and Robinson (1) investigated the 
vapor-liquid equilibrium of the nitrogen-butane system at 
temperatures from 100° to 300°r. ~ From “their "results; a 
retrograde region was found, and the K-factors for the 
nitrogen decreased rapidly as the temperature was increas- 
ed by increasing the temperature. They have drawn the 
critical locus over the range of data taken. They used 
the Edwards gas density balance to get the composition of 
their phases. Roberts and McKetta (53) reported data 
that considerably disagree with the results of Akers et al. 


They attributed the inaccuracy in the results of Akers and 
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co-workers to the extreme sensitivity required when ana- 

lyzing the phases with the density balance. Their study 

was conducted in the temperature range from 100° to 280°F 
and a pressure from 236 to 3402 psia. 

Poston and McKetta (42) experimentally determined 
the vapor-liquid equilibrium compositions for the nitrogen- 
n-hexane system over a pressure range from 250 to 5,000 
pSsia and at temperatures from 100° to 340°R. They found 
like other Jinvestigators (1,5, 53) that the solubility 
of nitrogen in the liquid hydrocarbon phase increases with 
pressure. 

Boomer, Johnson and Piercy (10) obtained the 
solubility and phase density data of nitrogen in n-heptane 
at 100.9 atmosphere and temperatures from 25° to 115°C. 
Akers, Kehn and Kilgore (3) reported the volumetric and 
phase behavior of the same system at temperatures between 
90° and 360°F. They found that at 90°F the density of the 
liquid phase increases with the increase of pressure. At 
higher temperatures, the greater solubility of the nitro- 
gen in heptane causes the density to decrease with increas- 
ing pressure. They attributed this unusual phenomenon to 


the non-ideality of the nitrogen-heptane solutions. 


Related Methane Binary Systems 
The methane-ethane system has been studied by some 
investigators (22). Price and Kobayashi (47) in their 


study of the vapor-liquid equilibrium of the methane- 
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ethane-pentane system reported their results of the 
binaries methane-ethane and methane-propane systems in 
the temperature range of -200°F to 50°F and the pres- 
sure from 100 to 1400 psia. No azeotropic mixtures were 
found in this’range. 

Sage, Lacey and Schaafsma (57) studied the methane- 
propane system and determined the densities of different 
mixtures in the pressure range of 10 to 200 atm. and tem- 
peratures from 20° to 90°C; the vapor-liquid equilibrium 
compositions were reported at pressures from 10 atm to 95 
atm, and temperatures from 20° to 90°C. They obtained a 
large region of retrograde condensation of the first 
kind. Reamer, Sage and Lacey (49) extended the range of 
study to a pressure up to 10,000 psia and temperatures 
From-240°"to 460°F) whtle Akers; Burns and Fairehiidr (2) 
investigated the methane-propane system at temperatures 
from -176° to 32°F and pressures from 50 to 1,450 psia. 
They drew the phase envelopes for various mixtures. 

Their data are in agreement with the results obtained by 
Sage, Lacey and Schaafsma. 

Neberbragt (38) in his study of the vapor-liquid 
equilibrium of the methane-butane system used n-butane 
containing a small percentage of iso-butane (96.6% n- 
butane and 3.4% iso-butane). He studied this system 
at pressures of 10 to 30 atm. and temperatures Erol ee: 


to 43°C. Sage, Hicks and Lacey (61) reported the 


= & - - 


- ‘edt Yo edivesx «told bostxoget moseye” ensineq-stedse 


mk amegaye easqorgq-smsrsom brs stsdso-sasdtom astusnid af 
-asig sid bas 1°02 oF TOOS~- to spmst sruts1agmet ots _— 
eiew Bexudxim oigorfosss of .sieq OOS oF OOL mozt owe 
-opnss eins nt Bayet - 
-6nardem oft beibuste (T2) smetssio® bas yeosd ,spse 
dneretIib to eeltienob odt benimzeteb bas moseye emsgozg 
tied bas «mgs OOS o¢ OL to spas stvaestg ont at asxudeio 
mrtzdifiieups brypii-xoqsy st .2°0R@ oF °OS mort asiuts18q 
2e@ oF mos Of mox? esrveestq ts betx1ogqser Sisw eroit LeoqgmoD 
& beatstde ysdT .5°0° o¢ °OS mort eexsdsisgmed bas \mos 
jaxi? so to dotd¢sensbnon sbs1rporte: to noips: sprst 
26 epas1 edt Bebaetxs (@)) ysosJ Bas spsa , 19msen pai 
gexrvssreqmss Bas steq 000,0f o+ qu siveeetdq 5 OF ybuse 
(S) DbifMotrsY bas enxwud ,2exsxA sltdw ,A°0ab oF OS mort 
@sivssreqnes +5 metaye sasqgoitg-ensnitem eds bstspiseevak 
-B6Leq 02@b,f oF ' 0e mor asrveesigq bas F°SE ot ° OV i=‘ mox : 
-Bsitvsxim agoitsy 10% asqolevns sesiq oft wath ysdT 
yd Bbenistdo adivuest sris dtiw tnemesips at sis S5teD <iodt 
-6metssdoe Bas yeond .9ps828 
bivpif-togsv sedi to ySuse aid ni (8£) tpsxdisden 
enstjud=-a beey mogeye onstud-sasdtem ont to mindi Livpe 
~m 3.30) sastud-o2f to spsinsoisg Iisme 56 et bi J 
meseye eins boibuste sH .(snstud-oet 8). bas onssad . 
“°LS= mOxt asistsxeqmet bas .m3s0€ oF OL to eorvesenq 35 


d ¥ 4 ' , : ety - ; 
9d3 betxoge: ([3) yoosd bos exci ,epse .d°ER e. " 
at) sy 7 


ie 


compositions and specific volumes of the cexisting phases 
of the methane-n butane system from 70° Vto 2507F and’ 31,3 
to 1923 psia. They showed the loci of the critical states 
and the cricondentherm. Roberts et al (54) extended the 
range of study from -80° to 280°F, their data are in 
agreement with Sage et al. Wang and McKetta (75) reported 
the vapor-liquid equilibrium data of the same system at 
-140°F and pressures from 73 to 453 psia. 

Sage, Webster and Lacey (62) studied the solubility 
of methane in n-pentane at 100°, 160° and 220°F at pres- 
sures’ up to 3000 psia. They found that the solubility of 
methane in pentane increases rapidly at high pressures. 
Boomer, Johnson and Piercy (8) determined the densities 
and compositions of the equilibrium phases of the n- 
pentane-methane system at 25% 55° and 85°C at total pres- 
sures from 35 to 190 atm. Similar measurements at 25°C 
and pressures from 35 to 135 atm. were made on a system 
composed of a mixture of n-pentane and iso-pentane with 
methane. The work of Sage and co-workers was found to 
agree with the work of Boomer et al., but a regular 
small discrepancy appeared which could be attributed to 
the presence of a small percentage of nitrogen in the gas 
used by Boomer and co-workers. Boomer et al found that 
methane is more soluble in the mixture of pentanes than 
in pentane alone; this is also in agreement with Hill and 


Macey. Taylor et al (73) reported the phase behavior of 
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the methane-n pentane system through the temperature range 
of 100° to 220°F, Their results indicate a progressive 
decrease in the critical pressure of this system with an 
increase in temperature, and a rapid change in the com- 
position of the dew point gas with pressure at the low 
pressure. Sage et al (60) determined the specific volume 
of six mixtures of methane and n-pentane for temperatures 
between 100° and 460°F, at pressures up to 5,000 psia. 
They reported the compositions of the coexisting phases 
at equilibrium as well as the equilibrium constants. 
They drew the loci of the cricondentherm, the critical 
state and the maximum two-phase pressure. 

Several investigators (41, 9, 19, 62, 66, 68, 70, 
24, 69) have studied the vapor-liquid equilibrium of the 
methane-n hexane system. Frolich et al (19) measured 
the solubility of methane in n-hexane at 25°C and at 
pressures as high as 1,300 psia, they found that the 
solubility of methane in normal hexane followed Henry's 
law fairly closely. Hill and Lacey (24) determined the 
solubility of methane in hexane at 86°F and 300 psia, 
while Sage, Webster and Lacey (62) went up to a pressure 
as high as 2,500 psia and temperature from 100" to 20S 
Their results, if interpolated to 25°C, will be lower 
than those reported by Frolich et al (19). Bubble point 
and dew point densities and COMPOST ELOnNS. ate 25.4, Joe sand 
85°C corresponding to pressures up to 230 atm., were 


reported by Boomer and Johnson (9) for the methane-n 
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hexane mixtures which contained a small amount of nitro- 
gen. The data reported by Boomer and co-worker is in 
Substantial agreement with that of Sage et al. Schooch, 
Hoffmann and Mayfield (66) determined the solubilities 

of methane in n-hexane at temperatures of 160° and 220°F 
and corresponding to pressures up to the critical. They 
showed that the maximum critical pressure exists at a 
temperature somewhere between 100° and 220°F. They also 
determined the densities and the compressibilities of the 
resultant solutions up to pressures of 6,000 psia. 
Stepanova and Vybornova (70) investigated the methane-n 
hexane system between 50 and 210 atm., and 0° and 60°C, 
they reported the critical pressure and equilibrium con- 
stants. Shim and Kohn (68) investigated the vapor-liquid 
equilibrium compositions of the methane-n hexane system 
at temperatures between -110° and 150°C and pressures up 
to 160 atm. They found that there is a gradual decrease 
in the solubility of methane in the liquid phase of the 
methane-light paraffin hydrocarbons, as the paraffin 
molecular weight increases. N-hexane is the highest 
molecular weight paraffin which is completely miscible 
with methane. The temperature range of 100° to 340°F and 
pressure range of 250 to 2,902 psia was covered by the 
study of Poston and McKetta (41). They showed that the 
solubility of methane in n-hexane decreases with tempera- 


ture. Their data are in agreement with the data of 
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Boomer (9), Sage (62), Schooch (66) and Shim (68) except 
for the critical temperatures reported by Shim and the 
vapor compositions reported by Boomer. Small disagree- 
ment with Boomer's composition data was also reported by 
Shim. Poston's data agreed to within 0.5 mole percent 
with data interpolated from Shim. 

Boomer, Johnson and Piercy (10) studied the methane- 
heptane system and reported the densities and compositions 
of both equilibrium phases at temperatures 25°, 55° and 
85°C, and pressures from 35 to 250 atm. Reamer et al (48) 
conducted their volumetric and phase behavior studies at 
pressures up to 10,000 psia in the temperature interval 
between 40° and 460°F. He found that at a temperature 
lower than -114.7°F, two liquid phases exist which is un- 
usual for paraffin hydrocarbon system. Kohn (29) made 
his study at pressures up to 1,500 psia and from -200° to 
S40 eRe The results of Reamer when compared with those 
of Boomer showed a large difference; while an excellent 
agreement was obtained between the data of Reamer and 
thatacft tKohn. 

Kohn and Bradish (30) determined the vapor-liquid 
compositions and molar volumes of the methane-n-octane sys- 
tem at eight different temperatures between -50° and 150°C. 
Temperatures, pressures and compositions along the three 
phase (vapor-liquid-solid) lines have been reported at 


temperatures between -107.6° and Sow. Cee lL ney, LOUund= tha & 


an a 
1) ; i. -_ _s my iu 


Lt a a oe sae 
ae . v ie nk: 


Ty LG? 


: teal mide bas (88) Mogorio® (Say opse. ( ) x 
end bas mide yd betiogez aarytersqmet sant 

ae ifsm2 .xemooa yd bodxogay enoitiaoqnon 30 

| ud bettoges oels 2sw SIBH obs teoqmos 2't9mo0c0d sie 


tneoteq elom 2.0 airntiw of beetes sisbh e'noseod «. mt 


.miie mort betslogzesnt sdeb itv 


~enaijem eit beibyta (Of) yoxsid das noeadol ,xesood “e Z 


— 


anoisteoqmoo bas selitiensb eid besroqet bas. meseys age [ 
bas °@@ \°GS sexwssregne: ts soasdq mutsdiltups died fo 


(8b) fs te temesh mts O28 oF CE moxk sexueaeuq Gas 44 ag 
; ¥ ne 
38. @eibuse toivaded easriq brs zirsemuiov ited Sesow 


fevteint exysexeqgnes eds nt saieg 000,0f o3 qu setues 


extuvsiteqomes © 3s gacd baveot sH ,.T°03S Bas °OR « 


= 


“mu ai doiriw Jaixp esesiq biupil ows ea°l  bLLH- ners te of 
ebem (€S) acon meteye sodissotbyd aliisisg 10} aes 
oF °O0S- mox® Buna aiag 002,1 ot qu setueesiq 35 ainda 
. ssonis ddiw bexrsqmoo serdw asmsef to asjlvesx oat at abe 


ae 
tnelfLeoxe as elinw ysonerst2ib episi s bswoda xemoof 20 
hues 


Bas temseh to s35h os nsewssd benissdo esw tneme Ips 


iad 


Adon tos 

Nilichdiinetiies ed3 benimxeteb (0€) sdetbsxd bas «dom 
“aye ensito-n-sassitom sdt to asmutov islom bas enok: 
span, bas °“O0¢- neewted eexutsreqmes tne1s?t ib stood 


; i 
~4 


| xy. a ial adt ehois eneatiecgae: bas a machcaae, a deem 


a oc bas °3. a neewded 
: ‘ “a ~s 


- 28 - 
Raoult's law does not hold accurately fon en—-octane= past 
0.20 mole fraction methane. 

Sage, Lavender and Lacey (59) studied the volu- 
metric behavior of several mixtures of methane-n-decane at 
pressures up to 4,500 psia and in the temperature inter- 
val between 70° and 250°F; the composition of the dew 
point gas was also determined at pressures up to 2,500 
peta at 100°,,°160°Safd®2202F .° "Reamer et®al-150)" in®his 
study covered the range of pressure from 20 to 5,300 psia 
and temperatures from 100° to 400°F. Koonce and Kobayashi 
(31) experimentally determined the solubility of methane 
in n-decane from measurements of pressure, temperature, 
and volume only. They calculated the K-values from methane 
in the methane-n-decane mixture at temperatures between 
-20° and +40°F, on the assumption that the vapor is free 
from the non-volatile components. Their results are con- 
Sistent with the ‘work of Sage and Lacey’ (58) carried out 
at higher temperatures. Beaudoin and Kohn (6) studied 
the vapor-ligquid equilibria of the methane-n-decane sys- 
tem at temperatures between -25° and 150°C; at pressures 
up to 100 atm. They found that Henry's law expressed as 
fugacity-mole fraction ratio holds as a reasonably accurate 


approximation at pressures below 30 atm. 


Related Ternary Systems 


Roberts and McKetta (51) experimentally determined 


the equilibrium compositions for the methane-n-butane- 
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nitrogen system between 100° and 250-8, anid 5 Uma ices) 000 
psia. They found that at 100°F the K-values of methane 
decrease with increasing the nitrogen concentration at 
500 psia, but this trend is reversed above 1,500 psia. 
This behavior was attributed to the difference in solu- 
bility of the methane and nitrogen in the liquid phase. 
As the temperature and pressue increase this difference 
becomes smaller and smaller, and the K-factors of methane 
increase with increasing concentrations of nitrogen at 
constant pressure. 

Boomer et al (9, 8, 10) studied the vapor-liquid 
equilibrium of the nitrogen-methane-pentane, nitrogen- 
methane-hexane, nitrogen-methane-heptane systemsat pres- 
sures from 35 to 250 atm. and in a temperature range from 
ZeetOmoo. Cor lney 1ound that the behavior of «the equa— 
librium constants is very similar for the three systems, 
with a regular variation from compound to compound. At 
all comparable pressures, the K-factor of the heavier 
hydrocarbon in the systems decreases when going from pen- 
tane to heptane, and changes most with temperature in the 
case of pentane. While the equilibrium constant for 
methane, at pressure above 50 atm. increases in going from 
pentane to heptane. The region of retrograde condensation 
is largest for heptane and least for pentane; and LSeLe= 
duced by increasing the temperature. 

Azarnoosh (4) investigated the methane-n-~decane- 


nitrogen system between 100° and 280°F, and 1,000 and 
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5,000 psia with the nitrogen concentration varying 

from 0% to 100% in the vapor phase. He found that the 
change in convergence pressure of this system is linear 
with the change of nitrogen concentration in the vapor 
phase since the mole fraction of n-decane in the vapor 
phase is extremely low. The convergence pressures at 
zero and one hundred percent nitrogen in the vapor phase 
are equal to the critical pressures of the binary systems 
of detane-methane and decane-nitrogen respectively at the 
same temperature. He found that as the mole fraction of 
nitrogen in the vapor phase increases, the K-factors of 
methane and nitrogen increase, while the K-factor of 


decane decreases, 


Multicomponent Systems 


Stutzman and Brown (72) experimentally determined 
the compositions of the vapor-liquid equilibrium phases 
for the common components of a natural gas including ni- 
trogen and methane, at a total pressure of 100 psia and 
temperature which varied from -130° to 240°F. Davis et 
al (15) presented the phase and volumetric properties of 
some natural gases and mixtures of natural gases with 
methane and nitrogen, trom —200°F to above the critical 
points. Jacoby and Rzasa (26) determined the vapor—-liguid 
equilibrium compositions for nitrogen, methane, carbon 


dioxide, ethane and hydrogen sulphide in a natural gas- 
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condensate system. Two different overal mixture compo- 
Sitions were investigated over the pressure range 500 
to 4/000epsiaratsl00°wal502 and) 200°R. Cosway and Katz 
(14) studied the systems containing hydrogen, nitrogen, 
methane and ethane at pressures of 500 and 1,000 psia 
and at temperatures of -100° and -200°F. They reported 
the compositions of the conjugate phases at equilibrium 
for three ternary and the quaternary systems. Lehigh 
and McKetta (32) reported the compositions of the conju- 
gate equilibrium phases for the nitrogen-ethane-butane 
system at. temperatures® from: 100°) to.280°F and pressures 
from 500 psia to the convergence pressure. 

Roberts and McKetta (52) studied the effect of 
Netccogen on, the vapor-liguid equilibrium ratios of hydro- 
carbons in natural gas condensate system. Their experi- 
mental work was carried out under constant temperature 
eonditions at 100°Fr and at 220°F, and at ditferent pres- 
sures between 500 and 3,000 psia. The nitrogen concen- 
tration in the vapor phase was varied from 10% to 50% at 
each pressure. They found that the K-factors of the 
hydrocarbons were lowered by the addition of nitrogen and 


the effect was more pronounced at 220°F than at 100°F. 


Nitrogen-methane-n hexane System 


Poston and McKetta (43) studied the phase equili- 
bria of the nitrogen-methane-n-~hexane system at tempera- 


tures from 100° to 340°F and pressures from 500 to 5,000 
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psia. They found that at 100°F and low pressures 

(around 500 psia) the addition of nitrogen to the system 
at constant pressure causes methane to be preferentially 
dissolved in the liquid phase. The liquid phase at low 
pressures is composed principally of n-hexane. Therefore, 
methane at low pressures is more soluble in n-hexane than 
is nitrogen. However, at pressures of 1,000 psia or more 
some Of the physical properties of methane and nitrogen 
approach one another and their solubilities in the liquid 
phase follow the same trend. Because of this the K-values 
of imethane “increase with “the “addition *of nitrogen in the 


vapor phase. 
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IV, -EXPERIMENTAL, STUDIES 


1. Experimental Apparatus 


(a) Equipment 


A very detailed and complete description of the 
Original design of the vapor-liquid equilibrium apparatus 
has been presented by Sandercock (63). The equilibrium 
cell used in these experiments was a 600 ml cylindrical, 
windowed cell of variable volume made of 316 stainless 
steel. It was designed to withstand a pressure of 8,000 
psia at a maximum temperature of 300°F. The freezing 
point of mercury, the confining liquid, determined the 
minimum operating temperature. The cell was designed to 
Contain»three sections: Sagas sections withsa floating 
piston, a restriction with a window assembly, and a liquid 
section having less capacity and smaller diameter than the 
gas section. 

the? positronlof the vapor=laguidianterféace: iny the 
cell could be changed without changing the volume of the 
equilibrium mixture.°-°This was accomplished» by connecting 
tive ‘top and: bottom Of thet Celie tostwormends. oftas double act— 
ing spump) andyain thismway. mercury icouddsbe;memoved si rom, one end 
at the same rate that it was added to the other end. fThus, 
the equilibrium mixture being confined between mercury and 
the movable piston, could be movedcr toiche desired level. 


The displacement pump connected to the bottom of the cell 
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was used to add or withdraw mercury thus decreasing or 


increasing the working volume of the cell. 


In order to see the inside of the cell around the 
sampling port, a periscope arrangement was used. Light 
Wasy>Drovideds throughs daglein.-by) Leined "buldsscyeretypeucs 
window and viewed through a telescope adjacent to the 
Opposite window. 

Mixing of the cell mixture was done by spraying 
the liquid phase and some mercury into the vapor phase 
using the double acting pump. A valve, placed in the low- 
er section of the cell, and having a small aperture pro- 
Widedpa Linewspray ofathe fluids) (65). 

Charging and sampling were done through a port at the 
window level. Samples were taken by means of a micrometer 
sample pump, having a total volume of about two cubic 
centimeters. In the original design, the sample pump was 
placed horizontally, but later (64) it was changed by a 
new pump placed vertically to prevent the possibility of 
a slug of the sample fluid being trapped in the pump. 

The cell and sample pump were immersed in separate 
ethylene glycol constant temperature baths. The pump 
was usually maintained at a higher temperature than the 
cell to avoid condensation of n-hexane during sampling. 

A bimetallic strip-controller was placed in an auxiliary 
bath, and used to govern the heater and refrigeration 


units giving a constant temperature for the ethylene 


elycolmcell bath. 
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A schematic arrangement of the experimental equipment is 
shown iin Bigure 1); 

The temperature of the cell was measured using an 
iron-constantan thermocouple inserted into the side of the 
cell. The thermocouple was calibrated in the temperature 
range from 32° to 200°F, against a standard precision 
mercury thermometer readable to 0.05°C, and a millivolt 
potentiometer, model 8686 manufactured by Leeds and Northrup 
Company, was used to measure the potential. The temperature 
was controlled to within +0.3°F and read to ‘within +0.1°F. 
The pressure of the system was measured using two Heise 
Bourdon gauges rated at 5,000 psia and 1,000 psia and 
marked at intervals of 10 psi and 2 psi respectively. Both 
gauges were calibrated with a Ruska dead weight tester. 

The pressure was adjusted manually within +0.5 psi up to 


1,000 pSiagand <torvwi thin +3 epsisaboveathat. 


(b) Leak Detection 

The equipment used for the study of the vapor-liquid 
equilibrium at high) pressure must be completely eiree Crom 
leaks. As a first step in the experimental program, there- 
fore, the equipment was pressure tested, and.all leaks were 
carefully eliminated. Leaking connections and valves were 
tightened or replaced. Studies including hydrogen sulphide 
had previously been conducted with the same equipment. 
This resulted in the deposition of mercuric sulphide on 
several parts including the valves, the pistons and cylin- 


ders of the pumps, and the wimlows, pistons and walls of 
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the cell. Several attempts were made to remove these 
deposits and it was finally found that sodium sulphide 
(39) was the most effective cleanser. 

In the course of the experimental program, this 
equilibrium equipment was used at working pressures as 
high as 2,000 psia. Although when pressure tested, the 
equipment could hold a pressure of 2,500 psia for several 
hours, many leaks appeared during the experimental runs. 
Most of the leaks in the cell were around the windows, due 
to flowing or breaking of the gaskets. It was concluded 
that the gaskets suggested by Sandercock (63) and success- 
fully used at pressures of about 1,500 psia, were not 
Suitable. Gaskets made of a harder material had to be 
used and after several trial and error procedures it was 
found that glass-filled teflon was a suitable material. 

Before starting the experiments, the sampling sys- 
tem was tested for leaks with a Veeco Mass Spectrometer 
Helium Leak Detector, Model MS-12AB-R. All leaks were 
eliminated. However, during the experimental work many 
leaks occurred in the sampling pump, most of them were at 
the packing, due to therflowing of the teflongpacking 
rings suggested previously (64). The pump operated for 
a much longer period with no leaks at the packing when 
gaskets made of teflon and similar to Chevron packing 


rings were used. 
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2. Materials Used 


The normal hexane used in the experiment was 
Philipps Petroleum Company research grade n-hexane with 
a purity of 99.99 mole per cent. The impurity determined 
by GLC analysis was methyl cyclopentane. Matheson of 
Canada supplied the research grade methane which had a 
purity of 99.97 mole per cent. 

Nitrogen was supplied by Canadian Liquid Air and 
had a minimum purity of 99.7 mole per cent. 

Clean triple distilled mercury was used. The pro- 
cedure used to clean the mercury has been detailed by 


Hughes (25), 


3 Experimental Method 


Charging the Cell 


The equilibrium apparatus was first evacuated over- 
night and then completely filled with mercury and pressure 
tested. To charge the equilibrium cell with the three 
components, hexane, having the lowest vapor pressure, was 
introduced first, then methane and finally nitrogen. 
Hexane was charged to the cell either by means of the 
sampling pump, or it was allowed to flow freely into the 
cell. In the latter case, some nitrogen was first intro- 
duced at the sampling port and some mercury was removed 
from the bottom of the cell, vacuum was then applied to 


evacuate the nitrogen, while the piston was at its lowest 
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position in the cell. Hexane was then admitted. The gases 
were charged to the cell by connecting the cylinder con- 
taining the gas to the sampling pump. The lines connecting 
the cylinder, pump and cell were evacuated; the gas was 
then allowed to enter the cell. When the pressure in the 
cell was equal to the pressure in the gas cylinder, the 
displacement pump could then be used to introduce more gas 


by removing some mercury from the bottom of the cell. 


Attainment of Equilibrium 
After the cell had been charged and contained the 


three components in two phases, the temperature and pres- 
sure were adjusted to the working conditions. The double 
acting pump was used to lower the liquid level below the 
internal valve. The double acting pump was then actuated 
rapidly enough to push the liquid through the internal 
valve, thus producing a fine spray into the vapor phase. 
This produced mixing of the phases and a change in pres- 
sure was noticed when the phases were not in equilibrium. 
The displacement pump was then used to adjust the pressure 
in the cell. The phases were allowed to separate and the 
mixing procedure was repeated until the pressure stabi- 
lized. The contents of the equilibrium cell were allowed 
to stand undisturbed for a period of two to four hours 

for the phases to completely separate. The liquid and then 


the vapor were sampled and analyzed. 
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Sampling Technigues 


The solution of the problems encountered in sampling 
the equilibrium fluids required considerable thought and 
experimentation. The aim of sampling is to transfer a 
small representative sample of the equilibrium phases 
from the cell to the chromatograph to be analyzed. 

The sample must represent accurately the composition 
of the phases in the equilibrium cell and they must be in 
the vapor phase while being injected into the chromatograph. 
The method used was to transfer a small sample at high pres- 
sure from the equilibrium cell to the sampling pump, then 
expand this sample into an evacuated container. A portion 
of the expanded sample was then injected into the chromato- 
graphi. 

(a) Sampling, Procedure 

Referring to Figure 1, the sampling procedure may 
be stated as follows: 

1. The sampling lines, pump and container were first 
evacuated with valves number 5, 21, 23 and 22 complete- 


ly open, and numbers 19, 20, 11 and 4 closed. 


2. The vacuum system was closed off at 5. 
3. Valve 4 was then cautiously opened and mercury allowed 
to fill the tubing between valves 4 and 19. Some mer- 


cury was allowed into the sampling pump while the 


phase to be analyzed was at the window level in the 


cell. 
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4. Valve 21 was then closed, and valve 19 opened while 
mercury was slowly pushed into the cell by the sampl- 
Ing pump. 

5. When the first droplet of mercury was observed through 
the window, the phase was ready to be sampled. Great 
care had to be exercised during sampling, in order not 
to disturb the equilibrium condition by changing the 
pressure of the mixture in the cell. Therefore while 
Slowly removing the sample with the sampling pump, 
mercury was added at the same rate from the bottom of 
the cell with the displacement pump, thus keeping the 
pressure in the cell constant. 

6. When enough sample had been taken, valve 19 was closed. 

7~ Ihe=vacuum"was*then cut off from: the sampling contain— 
er by closing 23. The sample was then expanded into 
the container by opening 20, and pushing the piston 
of the pump all the way up. Valve number 4 was then 
closed, and valves 11 and then 21 opened allowing 
helium to push mercury and any sample remaining in the 
lines into the container. Valve 22 was then closed. 
In the case of liquid samples, after taking the appro- 
priate volume with the sampling pump, the double act- 
ing pump was used to provide mercury at the window 
level, and a few milliliters of mercury were then 
taken out of the cell at the sampling port. 

As the samples were always diluted with helium gas 


before being injected into the chromatograph; a few trials 
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were made to determine a reasonable sample size and helium 
pressure, and it was found that good chromatographic res- 
ponses were obtained when the high pressure vapor sample 
was about 1.4 ml and the liquid about 0.3 ml; while the 
helium pressure was about 30 psia. 

To make sure that the liquid sample was vaporized 
before it was injected into the chromatograph, some calcu- 
lations were carried out to decide the minimum temperature 
to be used. Appendix A describes the calculations. The 
sampling lines and container were wrapped with insulated 
heating tape, while the sampling pump bath was heated 
using an immersion heater; to make sure that the sample 
temperature is always higher than 100°F. 

(b) Sampling Containers 

A few attempts were made to develop the sampling 
method and @cide won the sampling container. The sampl- 
ing container had to be of the proper size and shape to 
facilitate the transfer of the sample from the equilibrium 
cell to the chromatograph. Mention could briefly be made 
of some of the containers tried. 

The first trial was done with. a glass pear-shaped 
flask with replaceable silicon-sealed rubber at'*the top, 
two-way stopcock at the bottom and a three-way stopcock 
connected to the side of the flask. The sample was 
introduced mto the flask from the side stopcock (number 


22 in Figure 1), the helium from the bottom stopcock. 
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The sample was then taken by means of a syringe through 
the rubber stopper for subsequent analyses. The sample 
had to be injected directly into the column of the chroma- 
tograph by means of the syringe. An improvement of this 
arrangement was made by using another flask, having a neck 
at the top, with a replaceable silicon-sealed rubber 
stopper and a three-way stopcock at its side. Here the 
Sample was introduced first, followed by helium through 
the side stopcock. This flask had the advantage over the 
previous one that the helium introduced through the side 
stopcock washed any remaining sample and mercury in the 
sampling pump and lines into the flask. The samples were 
taken for analyses the same way as before. However, the 
glass sphere connections and stopcocks of these flasks 
could not stand the helium pressure without leaking; and 
also during injection there was a possibility of hexane 
condensing in the syringe. 

To overcome this difficulty, a steel bomb was used. 
Two stainless steel valves were connected to the bomb, 
one also connected to the sampling pump and the other to 
the vacuum line. After evacuating the bomb; the sample, 
mercury and helium were introduced into it. The bomb was 
then disconnected and well shaken for proper mixing. 
After proper mixing was achieved, the samples were taken 
for analyses through a hose connection which was attached 
to one of the valves. The injection of the samples was 


made into the loop of the chromatograph. 
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However, the possibility of hexane condensing in 
the injection lines was still present, and to avoid this 
a modification of the injection method was made, by 
connecting a hypodermic needle to one of the valves of 
the bomb allowing direct injection of the sample into the 
chromatographic column. This proved to be very efficient. 
The sampling container used during this work was a stain- 
less steel bomb of 500 ml capacity, and having a diameter 
Of aboute2 1/4 inches. Two 1/4 inch Hoke stainless steel 
valves (number 23 and 22 in Figure 1) were connected to 
the end of the bomb. 

(c) Injection Procedure 

Great care was taken during the injection of the 
sample from the bomb into the chromatograph. The injec- 
tion procedure can be summarized as follows: 

1. The sample was flashed into the bomb from the sampl- 
ing pump and the helium was admitted through valve 
number 11 in Figure 1 in order to push mercury and 
any sample out of the lines as mentioned above. 

2. The bomb was disconnected from the sampling arrange- 
ment and the wrapped heating tape was then removed. 

tea kiwvacuum was applied at the connection near valve 22 
tontiree the autoclave: coupling Of any droplet or 
mercury that might have been left after sampling. 

4. The bomb was then shaken, and good mixing of the con- 


tents of the bomb was assured due to the presence of 


mercury. 
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5. A hose connection, to which was soldered a fine hyper- 
dermic needle, was connected to valve 23. 

6. The bomb was then placed in a heating oven maintained 
at a temperature’d about 170°F: 

7. A few minutes later, the bomb was removed from the 
oven, valve 23 was opened, and the sample allowed to 
purge into the atmosphere for about one minute, this 
insured that no air would be left in the valve and 
hypodermic needle, At the same time, hot air was 
blown over valve number 23 and the injection needle, 
to prevent any condensation of hexane due to the 
threttlLing ‘effect. 

8. A sample was then allowed into the chromatographic 
column 1 as shown in Figure 2 through a silicon seal. 
ft was, found that an injection ime vot about 2 to 3s 
seconds was enough to obtain a reasonable size of 
sample. 

9. After injection, valve 23 was closed, the bomb shaken 
and then placed into the oven for the next injection. 

The time that passed during injection while the 
sampling bomb was outside the oven was never enough to cool 
the contents of the bomb; the heat capacity of the stain- 
less steel bomb being relatively so much higher than the 
contents. After the analysis was complete, the bomb was 
removed from the oven, its valves were opened and the mer- 
cury drained. The bomb was then connected to the sampling 


system, and helium allowed to flow through valve 11 Figure 1 
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for about 5 minutes to insure complete purging of the 
bomb and the system. The helium flow was then cut off and 
vacuum applied. When the bomb and pump were thoroughly 


evacuated, the system was ready for sampling. 


Analytical Techniques 
(a) Analytical Instrument 


A Kromo-tog Model K-2 chromatograph manufactured 
by Burrell was used for analyzing the ternary mixtures. 
Two U shaped glass columns could be in a constant tempera- 
ture air bath, its own temperature adjusting device, and 
its own constant temperature gas measuring and injecting 
systems. 

The signal generated by the thermal conductivity 
detectors was recorded on a Moseley 7100B strip chart 
recorder, equipped with a disc chart integrator Model 229 
manufactured for Hewlett-Packard by Disc Instruments 
Incorporated. 

A filter bed packed with molecular sieve having a 
mesh size of 45-60 was installed between the helium 
supply tank and the chromatograph. The filter guards 
against drifting base lines by removing the impurities 
from the carrier gas, improves chromatography, and leng- 
thens the life of the detector. A small 4-way Hamilton 
valve was installed on the chromatograph, to back purge 
by reversing the flow through the columns and thus extend- 


ing the life of the columns. 
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A switch was inserted, allowing one to reverse the 
polarity of the recorder; and thus making it possible to 
record any peak while the carrier gas was flowing in the 


reversed direction. 


(b) Chromatographic Methods 


Due to the wide difference in boiling points between 
nitrogen and n-hexane, it was impossible for one column to 
make separation of the three components studied. The columns 
needed had to give a good separation of the components in 
a short time. A few columns were tried and some are men- 
tioned below. 

Starting with two one meter columns, they were packed 
and conditioned. The first column was packed with 25% sili- 
con 550 on Celite, and the second with molecular sieve. 

The first column separated nitrogen and methane as one com- 
ponent and hexane as the other. The second column gave 
good separation for nitrogen and methane. The idea was to 
connect both columns in series, so that after getting the 
peaks of nitrogen and methane, the flow of helium and the 
polarity of the recorder were reversed so as to get a peak 
for hexane that was retained in the molecular sieve column. 
However it was found that a simpler way to apply this idea, 
was by packing the top 80 cms of a 2.5 meters column with 
25% Si 550 on Celite and the rest with molecular sieve. 
Applying the same principle as above, a very good separa- 


tion was obtained, but it was not possible to remove the 
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tailing obtained after reversing the flow. The peak 
corresponding to hexane was considered too wide and the 
length of the Si 550 portion insufficient, resulting in 
the tailing observed. 

Since a standard 2.5 meter glass column was used 
in the chromatograph, it became necessary to find a 
material which when packed in a reasonable short length 
would be capable of separating methane and nitrogen, and 
therefore giving enough length for Si 550 to separate 
hexane. Activated charcoal was found adequate for this 
Operation and the next trial was made with a column packed 
With ot 550 70n- Celite for the first 220.cems portion and 
about 30 cms of activated charcoal. This column gave a 
good separation of the components, but after about two 
months, it was found that the response of the chromato- 
graph had changed and no longer gave reproducible results 
on consecutive days. As some deactivation of the char- 
coal seemed to have occurred, it was found preferable to 
change the column. 

The next trial was to investigate the use of two 
Separate columns, as if they were in two diiterent chroma— 
tographs, and then to quantitatively relate the chromato- 
grams which were obtained from both columns. Making two 
injections for every analysis; injecting a sample ANTS 
the one meter column packed with Si 550 on Celite gave one 
peak for both methane and nitrogen and another peak for 
hexane. Then on injecting another sample into the other 


column which was packed with Si 550 on Celite for a length 


ms = i 7 _ 
a ns 
oi 


5 
een | sor ents Drlpxotro% miele: a4 
Schad oot betxeeitenoo aBW onsxod oF ‘entt .% 


mk pntdiuest ,ineloriived: noitxoq 022 12 ant aha pat’ : 
.bevieedo pail 


born esy nmuloa sesip todom ¢.° bushbnste 5 soaté 

& bri? O+ yissasDet amsosol Fi :igertpotsmoxtio abd: a] 

tt “dageal pers ‘efdencedey sp mi batosg tect doirw Ist da ; F 
San vaoporiin bis snsiitem paijsisqsa to sidsgso ed btwow, 
etsisqse ot oa2 £3 102 capnel dpvons pnivip srotesens * 

eins aoe. oasupebs BPatvot esw Dsoorsdo bet eviJDA _sasaed _ 
bexosg meses 6 djdiw sbhem asw iatitst gxon sda bas nolseteqo, " 
bas moitizeg em OSS terri oft 10t otiieD no O@e 18 As ew ii 


S s¥ep nimufon erdT .i[sonxsnia betsvitos to amo OF suods #1 


ity 
a 


7 


owd tuods xetisi sed ,einsnoaiion ait lo noiststsqeae Ras 


-osemotd> srt to Sanogest sit tsilt bavot Asw FE fp 
eilwest sidiouboxrgs: svsp tzspnoi on bos bopnreris Ben 
=1sdo oft to aoissvitoseb smm2 2h .eysbh ovitussenos | : pk 
og eldexeterq Bavot asw di ,bexiuo:0 sved ot bemeer es 
| tumuios sag 

ows to » sat etepssesyint ot 26w isixzs gxen sAT» 
“init gaoxodtib owt ai sisw yorit ti 25 , anmuioo aed 
Sweatanealies eda sisfor Wisvisssistasup of ners Bas ve 


So “aM By , Ae 
| “erm Loo fjod, mo1 bentstdo sxsw alee 


ym 2 6 esi tosiai veteyfene Yieva ok « 
29 m0 Ha ca ati besa fimufog: seen 
é i fle £ 
eet “eapossitn bas POR. af 
90 eft oat afqms: Melting priyootnt no ned? oa 
pnb = to2 eFtLeD no 022 be dtiw boxosq eew dokdw 


- 49 - 


of about 80 cms and molecular sieve for a length of about 
170 cms nitrogen was separated from methane and a separate 
peak was obtained for each of them. Meanwhile, the hexane 
was still retained in the column and was eluted by revers- 
ing the flow. The problem was to find a common basis to 
relate the peaks obtained from both columns. A reasonable 
method was to inject a known volume of sample to each 
column and therefore relate the obtained results from both 
columns of the chromatograph. It was tried to fill up the 
loops connected to both columns of the chromatograph with 
the vapor to be analyzed at atmospheric pressure. The 
relative volume of the loops would give the factor needed 
to correlate the responses obtained from both columns. 
Using this method for a binary mixture composed of methane 
and hexane at two different equilibrium conditions (100°F 
and 1,000 psia, andeloo0°r and 500 psia) the results obtain 
ed did not coincide with those reported in the literature. 
The wlative volumes of the loops were found not to be 
always the same. This fact along with the extreme preci- 
sion. required for analysis made the use of this method 
impracticable. 

The next method tried was by connecting two columns 
in series and getting the complete analysis by one injec- 
tion in less than 7 minutes. The chromatographic arrange- 
ment is shown in Figure 2. Column 1 is the one meter 
column, packed with 25% Silicon 550 on Celite, while column 


2 has a length of 2.5 meters, about the first 80 cms of the 
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left leg are packed with 25% Si 550 on Celite and the 


rest with molecular sieve. 


1 was placed between the two columns. 


The working side of detector 


A very good sepa- 


ration was obtained by this arrangement. 


The column conditions were 
Colum 1 

Column length 
Detecreor-Current 

Detector Bath Temperature 


Column Temperature 


Column 2 


Column Length 
DEBECtEGr CUuLrreire 
Detector Bath Temperature 


Column Temperature 


25% Si 550 on Celite 
1 meter 

200 mA 

LOO 7'C 


Sh Osea @ 


25%. S41) 550 on Celite, and 


Molecular Sieve 
2.5 Meters 
200 mA 
SOG 


BUR 


Helium, the carrier gas, had a flow of about 32.7 mL/Man. , 


measured at the sample exhaust side of detector 2. 


(Gan Analytical; Procedure 


For the analysis of a ternary nitrogen-methane-n- 


hexane mixture, the procedure used was the following: 


The one meter column was calibrated for pure n-hexane, 
and the 1.5 meters column for pure nitrogen and methane. 
The sample of the ternary mixture was injected directly 
into the first column, giving a peak for both methane 


and nitrogen, and a peak for hexane. Thus the hexane 
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content of the sample could be obtained from its 
Cali bpratron: 

3. As ititcamenout of) therfirst colunmetheisampliemas 
passed through the second column where hexane was 
retained by the molecular sieve and two peaks could 
be obtained for nitrogen and methane. From these 
peaks one may obtain the amounts of nitrogen and 
methane in the original sample. 

4. The hexane retained by the second column was then 
washed out by backflushing. 

A Hamilton one microliter syringe was used to 
inject different known amounts of n-hexane into the one 
meter column and get the calibration. 

Known gas mixtures of Helium-Methane and Helium- 
Nitrogen having a Helium concentration of above 80%, 
were used to calibrate the columns for methane and nitro- 
gen. A weighing technique was satisfactorily used to make 
the mixtures. Different amounts of methane and nitrogen 
were injected, using the glass loops of the chromatograph. 
The responses obtained from both columns were then record- 


ed. 


The calibration of the columns is given in Appendix 
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VARRESULTS 


The vapor liquid equilibrium behavior of the 
ternary nitrogen-methane-n-hexane system was studied 
at 40°F and 100°F. The experimental work was done 
ateypressunenor 500, 1).000- ands.2,000 psianfor tbothk 
temperatures. The analyses of the equilibrium phases 
are shown in Table l. 

These data are shown graphically in the phase 
diagrams figures 3 through 9. 

To smooth the experimental results, the data 
were plotted on a Y-X graph for each temperature studied. 
The best curve was drawn through the points for each com- 
ponent and the data were smoothed while keeping the re- 
lations tY = £X = 1.0 satisfied. 

The nitrogen concentration in the vapor phase was 
arbitrarily chosen to compare the results. Table 2 shows 
the equilibrium compositions and K-factors obtained from 
the smoothed curves at different concentrations of nitro- 


gen in the vapor. 
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TABLE 1 


Experimental Data 


Temp. |Press. Compositions, Mole Fraction 


oF PSIA Vapor Phase Liquid Phase 
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TABLE 2 


Smoothed Ternary Data 


Nitrogen Methane n-Hexane 


40°F, 500 PSIA 
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TABLE 2 (Continued) 


Nitrogen Methane n-Hexane 


100°F, 500 PSIA 
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VI DISCUSSION 


Most of the available correlations used to pre- 
dict the equilibrium relations for hydrocarbon systems 
exclude nitrogen. However, the equilibrium ratios 
could be calculated for systems composed of nitrogen and 
hydrocarbons by the Hadden, Grayson and Winn Nomograph 
(21, 78) and by the NGPSA Engineering Data Book (37), 


and by the Chueh-Prausnitz method (45). 


Comparison of K-factors with Winn and NGPSA Methods 


To calculate the K-factors by the Winn Nomograph 
andethe NGPSA, tthe critical Vocus, of vthernitrogen—n— 
hexane system was needed in order to estimate the con- 


vergence pressure of the ternary mixture. 


CrrEercal Loci 


To calculate the critical loci of the inzerogen— 
n-hexane and methane-n-hexane systems, the method pre- 
sented by Chueh (12) was used. 


The critical temperatures and the critical volumes 


of binary mixtures were given by 
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where 0 is the surface fraction, and is defined by 


275 
Xs Ser 
(6) = 
xX acy + x, ‘Tap 


where T19 and Vio are correlating parameters charac- 
teristic of the binary mixture. 

The critical pressure was calculated using the 
Redlich-Kwong equation of state with some alterations. 
The calculated critical temperature and volume were 
used in the equation, together with Ki5 which is a 
constant characteristic of the component l-component 2 
IncLeracLtion’. 

Three constants are thus needed to calculate the 


Critical locus Ci ga binary mixture... Table 3S shows “the 


constants used. 


TABLE 3 
See Jad) 
Binary Maxture = mer 7 7 12 
Nitrogen-Hexane -1.012 0.54 ORES 
Methane-Hexane Ons / OFZ5 0.08 


Poston and McKetta (42) gave the composition and 


the critical pressure of a nitrogen-n-hexane mixture at 


7 


Ap! Gee aN a9 , 
as « 


vi? : 
_- ¥ - 


ve onion ak are notjos13 soetwwe id ee “= 
ENS 


od rf 


aaa egtssemersq pnitsist109 ets ad bas aol 

. ,oxudxin yretid sd% to 
eit “an i. 
ue hie 
.enoissietis amos isiw etsta to noiisups paowl-r ADL. tt ny Ws 


edd poteyv fhetslusles esw oxveestq [spizizo 


e16w emuflov bas sxutsrogmet Lsoitixzs botsluoiso on? 
6 eL doidw oft dtiw xsdi2zspot \notseyps edz at 
§ Snendamos~I tnenogmoo oft to siteitssosaaddo 


ofa stsluolss oF bebsen eunt Sis atnstason ssait 


2 


@ao Swode € sidsT .sivsxim yisaid 5s to enool Lsoktito 


. beer ; 


€ aIaAT 
ot pe! I iS 24 Pig fF x = a 
7 + pe ae org tKiM 
=e 0 Ke a 12 cD” + 2. | 7 | 
8f.0 b2.0 Sto. L- ae 
80.0 @s.0 gayo~ we 


tae ’ 


o> sit avsp (sp) sa20H0m bas « 
ad~n-ngpoxsia s to suuaeesg 


- 65 - 


280°F. This experimental point was fitted and the con- 
stants shown in Table 3 were obtained. The constant used 
to calculate the critical temperature of the mixture was 
determined first from a knowledge of the critical tem- 
perature and composition of the experimental point. Ki 5 
was obtained by graphically extrapolating values obtained 
for other systems containing paraffin hydrocarbons and 
nitrogen. The parameter used to calculate the critical 
volume was then determined by means of the two other con- 
Stants and the experimental point. 

For the methane-n-hexane system the first two con- 
stants were obtained from curves available for paraffin- 
paraffin systems in terms of pure component critical 
constants, while Ki> was given by Chueh (12). 

Tables 4 and 5 give the results obtained, and 
Figure 9 shows the predicted loci, together with the ex- 
perimental points given by Poston. The composition of 
the critical points for the methane-n-hexane system was 
not experimentally available. 

From a knowledge of the critical locus of the 
nitrogen-n-hexane system, the convergence pressures 
of this binary were obtained at the experimental tempera- 
tures. Equation (17) was then applied to calculate the 
convergence pressure of the ternary system. The critical 


temperature of methane was used in this equation because 


the experimental temperatures were higher than the critical 
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The liquid phase composition obtained from the smoothed 
experimental data were used to calculate the convergence 
pressure of the ternary system. 

The K-factors were then read using the Winn 
Nomograph given by Hadden and Grayson. 

The equilibrium ratios for nitrogen, methane and 
n-hexane were also obtained at the experimental conditions 
by making some interpolation of values given in the NGPSA 
Engineering Data Book. The K-factors for each component 
were read from the book at the experimental conditions 
at convergence pressures of 2,000, 4,000 and 10,000 psia. 
Graphs were then constructed, where the K-factors were 
plotted against the convergence pressure using the pres- 
sure aS a parameter at different experimental temperatures. 
From these graphs the equilibrium ratios were obtained 
at the convergence pressures of the experiments. The 
experimental K-factors of the nitrogen-methane-n-hexane 
system are compared in Table 6 to the K-factors obtained 
from the Winn Nomograph and those from the NGPSA Engineer- 
ing Data Book. 

It could be said that the predictions made by the 
Winn Nomograph are better than those from the NGPSA. 


However, both methods cannot be accepted as giving good 
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predictions of the K-factors in this case. The reason 
could be that equation (17) should not be applied for 
this system, because the mass fraction of hexane was 
always high in the liquid phase, which made the change 
in the convergence pressure calculated by equation (17) 
quite small. Therefore, the effect of the change in 
the total composition was not well represented by the 
convergence, pressure. It gould then be concluded that 
this equation could not be applied to calculate the con- 
vergence pressure in this case, and that more data are 
needed to obtain a better estimate of the convergence 


pressure. 


TABLE 6 


Comparison of K-factors by Winn Nomograph and NGPSA 


K-factors with experimental K-factors 


Average Percent Deviation 
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Equilibrium Composition by the Chueh-Prausnitz Method 


To predict the bubble point temperature and the 
vapor composition of the nitrogen-methane-n-hexane 
system in the temperature range from 500°R to 800°R, 

a knowledge of the liquid composition and the pressure 
are’ needed. ~Prausnitz and Chueh (45) "describe the 
method” in detail.” The critical propertres’ anc other 
constants characteristic of the pure components are re- 
quired. For the possible binaries, these properties and 
constants are needed, together with the self interaction 
constants and the Henry's constants both given as func- 
tions of temperature. 

The binary parameters needed were either available 
or were determined during the calculations previously 
Made to predict the critical’ loci... The seli-interaction 
constant was not available for the binary nitrogen-n- 
hexane system. All the other relations were given by 
Prausnitz and Chueh. for a smaller range of temperature 
than the one of this work. They were extended to cover 
the experimental temperature range by the methods explained 
in the text. 

The HENRY'S program and the FITTING program 
given by Prausnitz and Chueh were used to determine the 
self-interaction constant for the nitrogen-n-hexane 
binary system as a function of temperature. The experi- 


mental data of Poston and McKetta (42) together with 
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those obtained in this work were used in making the cal- 
culations. Table 7 shows the results obtained. As shown 
in the table, the results predicted are not very accurate. 
This cannot lead to the conclusion that this method 
Should not be applied to predict the ternary equilibrium 
compositions, because some errors could have been intro- 
duced during the extrapolation of the Henry's and the 
self-interaction constants. “Aliso, the binary “experimental 
data used to calculate the self-interaction constant for 
the nitrogen-n-hexane system were not complete enough to 
give a truly representative value. 

However, from Figure 9 it may be concluded that 
the Chueh-Prausnitz method gives a good prediction of 
the critical locus for the binary systems containing 
methane and hexane. 

Appendix C contains a comparison of the temperature 
and the vapor composition of the nitrogen-methane-n-hexane 
system calculated by the Chueh-Prausnitz method with the 
experimental data given by Boomer and Johnson (8) and 
Poston and McKetta (43). A comparison of the results by 
Poston and McKetta for the nitrogen-n-hexane and methane-n- 
hexane systems are also included. 

It should be mentioned that while calculating the 
self-interaction constants oa for the nitrogen-n-hexane 
system by the FITTING program, the Corresponding dilation 


constants n were given in the output of the program. 
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Some values of yn and the values of a calculated from them 
were neglected before plotting and fitting this constant. 
The reason is that the model presented by Prausnitz and 
Chueh does not allow for negative: 7. 

The FITTING program gives intermediate values 
from which a could be calculated once n is known. This 
procedure gave values for a that were far from those 
listed for similar systems and therefore these results 


were not used. 
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NOMENCLATURE 

£ fugacity 

£° reference fugacity 

ie denotes a function of variables 

K equilibrium ratio 
Ky > characteristic constant for 1-2 interaction 
M compositon parameter based on mass fraction 
m mass fraction 
im total pressure of system 
p partial pressure of a component 
p? vapor pressure of a component 
Pa critical pressure 

Pom critical pressure of a mixture 
Px convergence pressure 

R gas constant 

- temperature of system 

To critical temperature of a component 

Tom’ Tor critical temperature of a mixture 

Vv molar volume of the vapor or liquid phase 
Vo critical volume 

om! Vcr critical volume of a mixture 

X mole fraction in liquid phase 

Y mole fraction in vapor phase 


Z compressibility factor 
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Greek letters 


interaction constant 

dilation constant 

surface fraction 

correlating parameter for critical volume 
correlating parameter for critical temperature 
activity coefficient 

solubility parameter 
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Appendix A 


To make sure that n-hexane in the liquid samples 
would be transferred to the sampling bomb in the gaseous 
state, and no condensation would take place, the following 


calculation was made: 
pV = NZRT 


where V = Volume of the sampling container, taken as 500 ml 
R = Gas constant, taken as 670 ml.psia/gm mole -e 
T = Temperature of the sampling container, OR 
p = Partial pressure of n-hexane, calculated in psia 
at Temperature T, assuming ideal gas behavior 
Z = Compressibility factor, taken as 1.0 at the 


pressure and temperature in the sampling con- 


tainer 
N = Moles of n-hexane sampled 
Specific gravity of n-hexane = 0.66 
Molecular weight of n-hexane = 86.2 


Assuming the volume of n-hexane taken from the equilibrium 


cell = 0.4 mlz 


NZRT 
V 


, (0.4) (0r66) (clea®)). 67.0) 29. 
86.12 


(on) 


0 


By calculating the partial pressure at different Veenpera- 
tures and from a knowledge of the vapor pressure of 


n-hexane, Table Al and Figure Al were constructed. 
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FIG. -Al- DETERMINATION OF THE MINIMUM ALLOW- 
ABLE TEMPERATURE OF THE SAMPLING 
SYSTEM. 
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O. 
P is the vapor pressure of normal hexane. 


Therefore if t, the sample temperature is always 
maintained higher than 100°F, liquid hexane will definitely 


be in the vapor state. 
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Appendix B 


Calibration of the Chromatographic Columns 


Column 1 was calibrated for n-hexane, while column 
2 was for nitrogen and methane. When the component to 
be analyzed was eluted from the column, it passed through 
the thermal conductivity detector and a peak was subse- 
quently drawn by the recorder. The area under the peak 
was measured by the disk integrator. The number of moles 
of material entering the detector cell was proportional 
to the area under the peak. The relation between the 
number of moles of material that was added to each column 
and the area under the corresponding peak was represented 
by a straight line down to an experimentally obtained 
limit. From: this limit to the origin, the function was 
represented by a quadratic equation which passed through 
the origin and intersected the linear relationship tan- 
gentially. Therefore, if Y = number of moles and 
X = area under the peak, the relation Y = CX + px? will 
hold for X smaller than the experimentally obtained limit, 
above which the relation Y = AX + B applies. 

Where C, D, A and B are constants for the component 
to be analyzed. 

The experimentally obtained calibration points were 


fitted using a least square method. The results are pre- 


sented in Table Bl. 
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Appendix C 
TABLE Cl 


Comparison of the Temperature and Vapor Composition Cal- 
culated by Chueh-Prausnitz Method with the Experimental 
Values for the Nitrogen-n-Hexane Binary (42) 
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